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Abstract 

Vitamin D is well known for its regulatory role in calcium and phosphate homeostasis, but its role in muscle mass and 
strength during growth remains inconclusive. We explored the association of serum 25-hydroxyvitamin D (25(OH)D) with 
muscle development in girls from 11 to 18-years old. Whole body lean tissue mass (LM WB ), appendicular lean mass (aLM), 
muscle cross-sectional area at the lower leg (mCSA), maximal voluntary contraction of elbow flexors (MVC e | bow ) and knee 
extensors (MVC knee ) were assessed in 217 girls aged 10-13 years (at baseline), 215 in 2-year and 226 in 7.5-year follow-up. 
Serum concentration of 25(OH)D and intact parathyroid hormone (PTH) were analyzed retrospectively and girls were 
categorized according to their 25(OH)D levels (consistently insufficient 25(OH)D G LL <50 nmol/l and consistently sufficient 
G H h >50 nmol/l from baseline to 7-year follow-up). We found that 25(OH)D level declined until menarche (p<0.05) while 
LM WB , aLM, mCSA, MVC e | bow and MVC knee continued to increase (p<0.001 for all) post menarche. At pre-menarche, the G LL 
(n = 34) had higher LM W b and aLM than the G HH (n = 21 , p<0.05), while post-menarche the G HH (n = 1 5) had a greater catch- 
up gain in LM WB (p = 0.004), aLM (p = 0.001) and mCSA (p = 0.027) compared to the G LL (n = 65) over the first 2-year period. 
At the age of 18, no differences in muscle mass/strength between the low (n = 151) and high (n = 77) levels of 25(OH)D 
groups were found. This finding was independent of vitamin D receptor genotype and other confounders. In conclusion, 
our results showed that levels of 25(OH)D have no significant negative influence on the development of muscle mass and 
strength during pubertal growth both with longitudinal and cross-sectional comparison. On the contrary, our results 
suggest that the temporary negative association between 25(OH)D and muscle mass arises as a consequence of fast growth 
prior to menarche, and this negative association is diminished through catch-up growth after menarche. 
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Introduction 

Recent epidemiological studies have shown that people with 
poor early growth have an increased risk of sarcopenia in later life, 
indicating that sarcopenia may have its origin during early 
development [1,2]. Pubertal growth with optimal exercise and 
nutrition is considered to be a prime time for maximizing one's 
genetically predetermined peak bone mass [3], but it is not clear if 
that is the case with muscle mass. It has been suggested that 
vitamin D has beneficial effects on the development and 
maintenance of muscle mass [4,5]. The beneficial effects of serum 
25 hydroxyvitamin D (25(OH)D) (an index of vitamin D status) are 
mediated by its receptor in skeletal muscle [6,7], and genotypic 
variations for this receptor have been reported to be associated 
with differences in muscle strength [8,9] . However, recent studies 
on the effects of 25(OH)D have primarily focused on elderly 
individuals who are most prone to sarcopenia [10,1 1,12]. Only a 
few studies have considered adolescents [13,14]. These previous 



studies are cross-sectional in design and thus unable to examine 
whether the relationship between serum 25(OH)D level and 
muscle mass/ strength is independent of growth and maturation 
[8-11,13-15]. 

Puberty is a critical period of growth and development that is 
highly dependent on adequate nutrition. Vitamin D is well known 
for its regulatory role in calcium and phosphate homeostasis, and 
deficiency in vitamin D in childhood can cause reduced bone 
mineral density, delayed growth, skeletal deformities and impaired 
motor development [16-23]. However, the role of vitamin D in 
muscle mass and strength during growth is still inconclusive. In 
this study, we used a 7.5-year longitudinal dataset of pubertal girls 
to test our hypothesis that 1) a consistently low level of 25(OH)D 
has a significant negative influence on the development of muscle 
mass and strength in girls during the rapid pubertal growth period, 
and 2) the associations between level of 25(OH)D and muscle mass 
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and strength gain are dependent on maturation stage and vitamin 
D genotypes. 

Subjects and Method 

Study Design 

The subjects and the study design have been described 
previously [24,25]. Briefly, 396 girls aged 10-13 years (mean 
age 11.2 years, at baseline) were recruited from local schools in 
the city of Jyvaskyla and its surroundings in Central Finland to 
participate in a longitudinal study of determinants of body 
composition during puberty growth (the Calex-study). During 
the first 2-years, girls participated in a randomized intervention 
with calcium, vitamin D and dairy products supplementation. 
After the first 2-years intervention, we continued to follow the 
subjects on average for additional 5 years. The same data 
collections were performed at the 4- and 7.5-year time points. 
Because there were no intervention effects on muscle mass and 
strength the data were pooled together for the secondary 
analyses [24,25]. 

Study Population 

Among the eligible participants who had both 25(OH)D and 
body composition assessments, 2 1 7 of them were included in the 
baseline, 215 in 2-year, 110 in 4-year and 226 in the 7.5-year 
follow-up (of those 226 girls, 101 had both baseline and 7.5-year 
follow-up assessments). At the 4-year follow, only part of the 
participants was randomly invited for body composition assess- 
ments due to limited funding. In other time point's data 
collections, the main reasons for dropout were relocation, loss of 
interest and lack of time. 

Background Information 

Current health status, chronic conditions and regular medica- 
tions were checked by study nurse and physician. Dietary 
information was obtained from a food-intake diary kept for three 
days (two ordinary school days and one weekend day) as described 
elsewhere [26] . For this current report, the nutrient data related to 
body composition was used. Leisure time physical activity was 
assessed by questionnaire in terms of times/week and hours/week 
of participation in exercise. The age at menarche was defined as 
the first onset of menstrual bleeding and was determined by 
questionnaire or phone call during the follow-up. Body weight and 
height were measured with subjects wearing light clothes and on 
bare feet. 

25-Hydroxyvitamin D (25-OHD) and Parathyroid 
Hormone (PTH) Assessments 

The laboratory tests performed within two weeks' period during 
the same month of the year throughout the 7.5-year follow-up to 
avoid seasonal effects. Blood samples were taken in the morning 
between 7:00 and 9:00 after an overnight fast. Serum samples 
were stored at -70°C until analyzed. Serum 25(OH)D concentra- 
tions were measured by radioimmunoassay (Incstar Corporation, 
Stillwater, MN). The intra- and inter-assay CVs were 10% and 
15%, respectively. The reference range for 25(OH)D was 25- 
120 nmol/L. Serum intact parathyroid hormone (PTH) concen- 
trations were measured by using an immunoradiometric method 
(Nichols Institute, Juan San Capistrano, CA), with 10-65 pg/L as 
the reference range. The intra- and inter-assay CVs were 4% and 
3%, respectively. 



Vitamin D Receptor Gene Polymorphisms 

Genomic DNA was extracted and purified from EDTA blood 
samples using the QIAmp Blood Kit (Qiagen GmbH, Hilden 
Germany). TaqI polymorphism at exon 9 of the VDR gene, Apal 
polymorphism at intron 8 of the VDR gene and a C to T variation 
located at the restriction enzyme Fokl site at the translation start 
site of VDR gene were analyzed [27]. Capital (T, A, F) and small 
(t, a, f) letters refer to the absence and presence of the restriction 
site, respectively. Genotypes AA, TT, and FF represent homozy- 
gotes of frequent alleles; Aa, Tt, and Ff represent heterozygotes, 
and aa, tt, and ff represent homozygotes of infrequent alleles. 

Body Composition, Muscle Mass and Size Assessments 

Body composition was measured by dual-energy X-ray absorp- 
tiometry (DXA, Prodigy; GE Lunar Corp., Madison, WI USA). 
The CV of two repeated measurements on the same day was on 
average 1.0% for whole body lean mass (LMwb) and 2.2% for fat 
mass (FM). Appendicular lean mass (aLM) was calculated as the 
sum of lean mass in arms and leg as the surrogate for muscle mass. 

In addition, the left lower leg was scanned using peripheral 
quantitative computerized tomography (XCT 2000, Stratec 
Medizintechnik, Pforzheim, Germany). The scan location was at 
60% of the lower leg length up from the lateral malleolus of the 
fibula. The in-plane pixel size was 0.59 mm x 0.59 mm. The 
muscle cross-sectional area (mCSA) was analyzed using validated 
software (Geanie 2.1, Geanie®CE, Commit Ltd, Espoo, Finland). 
A contour was drawn manually along the outer boundary of 
muscle to eliminate the subcutaneous adipose tissue before 
analysis. The threshold for muscle was 10-279 mg/cm 3 . The 
CV of two repeated measurements on the same subject on the 
same day was on average 1% for mCSA. 

Muscle Strength and Performance Assessments 

Maximal isometric voluntary contraction force (MVC, N) of the 
left elbow flexors and left leg extensors was measured in a sitting 
position with an adjustable dynamometer chair (Good Strength, 
Metitur, Jyvaskyla, Finland). Subjects were encouraged to exert 
their maximal effort during an isometric test for ~3 s. The MVC 
of the elbow flexors was measured with the subject's forearm 
strapped at the wrist using velcro straps with elbow at 90° and the 
thumb in an upward position. Seat height was adjusted to ensure 
that the upper arm was at a 90° to the chest and the elbow was at 
the pivotal point of the apparatus. In the measurement of MVC of 
knee extensors, the subject's lower leg was strapped at the distal 
end of the lower leg using Velcro straps with knee at 120°. Seat 
height was adjusted to ensure that the hip was at 90° to the trunk 
and the knee was at the pivot point of the apparatus. The hip was 
also strapped by Velcro straps to limit its movement. The CV was 
3.7% and 6.6% for elbow flexors and knee extensors, respectively. 

Ethics 

The study procedures followed were in accordance with the 
Helsinki Declaration of 1975 as revised in 1983. The study 
protocol was approved by the ethical committee of the University 
of Jyvaskyla, the Central Hospital of Central Finland and the 
Finnish National Agency of Medicines (memo 1/2005, 
19.01.2005, 3/2006, 5/4/2006, 20/2/2007, 18/12/2007, 22/4/ 
2008, 6/6/2008, and 22/8/2008). The participants provided their 
written consent in accordance with the guidelines laid down by the 
ethical committees. When the girls were under the age of 18 years 
old, both participants and their parents signed the informed 
consent prior to the assessments. 
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Figure 1. The growth curves of 25(OH)D (a), lean mass(£= LM WB & c=aLM), muscle size (</=mCSA) and muscle strength 
(e=WIVC e | bow & f= WIVC knee ). TRM = time (months) relative to menarche; LM WB = whole body lean mass; aLM = appendicular lean mass; 
mCSA = muscle cross-sectional area; MVC e | b ow = maximum voluntary muscle contraction of elbow flexors; MVC knee = maximum voluntary muscle 
contraction of knee extensors. Grey dots and lines indicate individual values and the black solid line indicates the best fitting of growth pattern by 
hierarchical modelling. 
doi:10.1371/journal.pone.0082124.g001 



Statistical Analyses 

All data were checked for normality using Shapiro-Wilk test in 
IBM SPSS 20.0 for Windows. A hierarchical (multilevel) model 
with random effects (ML wiN 2.22 software, Institute of 
Education, University of London, UK) was used to explore the 
patterns of longitudinal changes of 25(OH)D and muscle mass and 
strength during growth. The hierarchical model allows inclusion of 
the data from every subject regardless of irregularity of temporally 
spaced follow-up or missing data. Time relative to menarche 
(TRM), instead of age, was entered as the explanatory variable in 
the form of polynomial functions to explain the evolution of target 
variables. 

The associations between longitudinal changes of muscle 
parameters and changes of 25(OH)D pre- and post-menarche 
were assessed by the following hierarchical model: 

y-,j = fi0jjcons+ /? 1 x Ln25(OH)Dij x dumiA 
+ P2 x Ln25(OH)Dij x dumiJ. 

with fSOjj = j80 + uOj + eOy 

where the outcome variable y tJ is LM WB , aLM, MVC elbow and 
MVC km . e . The subscript indices represent the zth measurement 
occasion for thejth individual, respectively. The predictor variable 
is Ln(25(OH)D). P 0 is the intercept, and both u 0 j and e oi j form the 
random portion of the model, whose means are equal to zero. The 
terms dumi_l and dumi_2 are the dummy variables introduced to 
denote the pre- and post-menarche periods respectively. Here, the 
time of menarche itself is selected as a shift knot for the model, 
which means that the coefficients could be different on either side 
of this time point. Thus, the associations between 25(OH)D and 



muscle mass and strength can be assessed by regression coefficients 
(31 for pre-menarche, and by (52 for post-menarche, respectively. 

In order to investigate whether the level of 25(OH)D had a 
significant influence on the development of muscle mass and 
strength during pubertal growth, we categorized girls according to 
their 25(OH)D levels and compared the two groups: 
Gll = 25(OH)D consistendy lower than 50 nmol/1 both at 
baseline and 7.5-year follow-up, and G H h = 25(OH)D consistendy 
higher than 50 nmol/1 both at baseline and 7.5-year follow-up. T- 
test was used to compare the two 25(OH)D groups. 

In order to evaluate whether maturation status has influence on 
the relationship between 25(OH)D and development of muscle 
mass/strength, we used a mixed model to compare the percentage 
changes between G LL and G HH and between pre- and post- 
menarcheal girls controlled for age at the 2-year follow-up, and the 
percentage changes of body height, PTH, intakes of vitamin D and 
physical activity level. If there was overall significance (i.e. 
interaction between vitamin D group and maturation stage by 
measurement time points), pair-wise comparisons between vitamin 
D groups were assessed by the Sidak method within pre- and post- 
menarche. The results are reported as estimated means ± 
standard error (SEs). A P value of <0.05 was considered 
statistically significant. 

Analysis of variance (ANOVA) was used to compare the group 
differences in VDR genotypes. In addition, when comparing the 
group differences in muscle mass and strength VDR Apal alleles 
were controlled for and P-values for the pairwise comparisons 
within pre- and post-menarche analyses were made using the 
Sidak method. 

Results 

At the age of 1 1, 13 and 18 years, the proportion of subjects 
with vitamin D insufficiency [(25(OH)D >25 nmol/1 and 
<50 nmol/1] was 56.5% and 59.3%, 63.4% and vitamin D 
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deficiency [(25(OH)D <25 nmol/1] 6.0%, 12.6%, and 8.9%, 
respectively. 

The hierarchical model with random effects for longitudinal 
changes of serum 25(OH)D, LM WB , aLM, mCSA, MVCaw and 
MVCtnee in girls are presented in Figure 1 . The concentration of 
25(OH)D declined till menarche (p<0.05) but not post-menarche. 
LM WB , aLM, mCSA, MVC elbow and MVC^ continue to 
increase after menarche (p<0.001 for all). 

We found that during pre-menarche 25(OH)D was negatively 
correlated with LM WB (P = - 2.048 (SE = 0.477), p<0.05), aLM 
0= -1.468 (SE = 0.388), p<0.01), mCSA (P = -0.323 
(SE = 0.097, p<0.01 ), MVCei W (P= -0.041 (SE = 0.019), 
p<0.05), and MVC^ (P= -0.050 (SE = 0.026), p>0.05), while 
at post-menarche, these negative associations disappeared (LM WB 
0 = 0.196 (SE = 0.484), p>0.05), aLM (P = 0.705 (SE = 0.394), 
p>0.05), mCSA (P = 0.047 (SE = 0.099), p>0.05) MVC dbow 
(P = 0.030 (SE = 0.020), p>0.05), MVC knee (P = 0.025 
(SE = 0.026), p>0.05). 

Menarche is a critical time point for the development of girls. In 
order to evaluate whether maturation status has influence on the 
relationship between 25(OH)D and muscle mass/strength, we 
divided girls into pre- and post-menarcheal status at the first 2-year 
follow-up and compared the Gll and G H h groups. We found that 
the G LL was more mature as reflected by TRM (p = 0.033), and groups. The G LL tended to differ from G HH in Apal (AA/Aa/ 
taller (p = 0.034), and had lower PTH level compared to the Ghh aa G LL = 26/64/10% vs. G HH = 55/27/ 18%, p = 0.09). Control- 
(p = 0.0 28, Table 1). No other differences in terms of level of ling for VDR Apal, the G LL group had higher LM WB (p = 0.02), 
physical activity and dietary intakes of micronutrients between the aLM (P = 0.021) than G HH while no significant differences were 
G LL and G HH groups were found. Further, there were no found in mCSA, MVC elbow and MVC^ in pre-menarcheal girls 
significant differences in terms of distribution of different VDR but these differences were not present in post-menarcheal girls 
genotypes of TaqI and Fokl between the G LL (TT/Tt/tt = 4 1 /49/ (Table 2). 

10% and FF/Ff/ff= 28/64/8%, respectively) and the G H h (TT/ In addition, no significant differences were found between the 

Tt/tt = 42/50/8% and FF/Ff/ff= 17/75/8%, respectively) G LL and G HH in any of the variables except for 25(OH)D level 



Table 1. Comparison of anthropometry, physical activity, intake of micronutrients, and PTH among girls grouped according to 
whether their 25-OHD levels were consistently sufficient, insufficient, increased or decreased over the first 2-years follow-up, 
according to their menarcheal status (ANOVA; mean and SE are given, adjusted for multiple comparisons by the Sidak method). 





Variables 


At the 2-year follow-up pre-menarche 


at the 2-year follow 


up post-menarche 




G LL (n = 34) 


G HH In = 21) 


G LL (n = 65) 


G HH (n = 15) 


Age (years) 


12.6 (0.1) 


12.8 (0.1) 


13.6 (0.1) 


1 3.0 (0.2) 


TRM (months) 


-8.7 (1.4) 


-11.8 (1.9) 


10.1 (1.0) 


6.2 (1.9)* 


Height (cm) 


154.9 (1.0) 


150.7 (1.3)* 


160.3 (0.8) 


158.8 (1.3) 


Weight (kg) 


46.9 (1.7) 


42.8 (2.2) 


53.4 (1.3) 


52.4 (2.6) 


BMI (kg/m2) 


1 9.4 (0.6) 


18.7 (0.8) 


20.8 (0.4) 


20.8 (0.9) 


FM (kg) 


12.6 (1.2) 


11.5 (1.5) 


14.7 (0.9) 


14.3 (1.8) 


Exercise(times/wk) 


3.00 (0.4) 


3.54 (0.5) 


3.06 (0.3) 


3.38 (0.5) 


Exercise (hrs/wk) 


2.66 (0.4) 


3.15 (0.5) 


3.02 (0.3) 


3.73 (0.6) 


Intakes of Calcium (mg/d) 


969 (81) 


1068 (100) 


1041 (59) 


1057 (118) 


Phosphorus (mg/d) 


1205 (66) 


1243 (81) 


1210 (48) 


1270 (96) 


Vitamin D (jig/d) 


3.01 (0.4) 


4.15 (0.5) 


2.0 (0.3) 


4.37 (0.6)** 


Energy (kcal/d) 


1586 (176) 


1769 (136) 


1899 (90) 


1567 (176) 


Protein (E%) 


15.2 (0.6) 


15.6 (0.8) 


15.3 (0.5) 


15.5 (0.9) 


Fat (E%) 


31.6 (1.0) 


34.4 (1.3) 


33.2 (0.7) 


30.6 (1.5) 


Carbohydrate (E%) 


53.2 (1.3) 


49.9 (1.6) 


51.5 (0.9) 


53.9 (1.9) 


Serum PTH(pg/ml) 


32.1 (1.9) 


27.4 (2.5) 


30.3 (1.4) 


24.2 (2.9)* 



Pre- and post-menarche was divided by TRM value of the 2-year-follow-up time point where menarche is 0, — TRM is pre- and +TRM is post-menarche. TRM = Time 
(months) relative to menarche; BMI = weight/height 2 ; FM = fat mass. 
Significant difference between the G LL and G H h: *<0.05; **<0.01. 
doi:10.1371/journal.pone.0082124.t001 



Table 2. Comparison of lean mass and muscle strength 
between girls whose 25-OHD levels were consistently 
sufficient and consistently insufficient, according the 
menarcheal status at 2-year follow-up controlled for VDR Apal 
(ANOVA; mean and SE are given, adjusted for multiple 
comparison by the Sidak method). 







Pre-menarche 


Post-menarche 




G LL (n = 34) 


G HH (n = 21) 


Gll (n = 65) 


G HH In = 15) 


LM WB (kg) 


33.0 (0.69) 


29.4 (1 .06)* 


36.3 (0.5) 


35.8 (1.05) 


aLM (kg) 


14.2 (0.35) 


12.8 (0.54)* 


15.7 (0.25) 


15.7 (0.53) 


mCSA (cm 2 ) 


48.4 (1 .64) 


44.9 (2.54) 


55.4 (1.20) 


56.2 (2.51) 


MVC elbow (N) 


141.2 (4.5) 


135.8 (7.0) 


154.6 (3.3) 


162.0 (6.9) 


MVC knee (N) 


351.8 (14.2) 


353.3 (22.0) 


399.4 (10.3) 


432.4 (21.7) 



LM WB = whole body lean mass; aLM = appendicular lean mass; mCSA = muscle 
cross-sectional area; MVC e | bow = maximum strength of elbow flexors; 
MVC knee = maximum strength of knee extensors. 
*<0.05 pairwise comparisons between G LL and G H h- 
doi:1 0.1 371 /journal.pone.00821 24.t002 
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Table 3. Comparison of muscle mass and strength among 
girls grouped according to whether their 25-OHD levels were 
sufficient and insufficient at the 7.5-years follow-up (T-test; 
mean and SD are given). 





Variables 


G L „„ (n = 151) 


G H ig„ (n = 77) 


LM WB (kg) 


38.2 (4.1) 


37.8 (3.9) 


al_M (kg) 


17.0 (2.3) 


16.5 (2.1) 


mCSA (mm2) 


6418 (1085) 


6212 (1028) 


MVC elbow (N) 


165 (30) 


169 (40) 


MVC knee (N) 


388 (93) 


386 (81) 


doi:10.1371/journal.pone.0082124.t003 



when girls at the age of 18 years. Further, including whole cross- 
sectional sample who participated in 7.5-year follow-up tests 
(n = 226), the results remained the same (Table 3). 

We further compared the percentage change of the muscle mass 
and strength between the Gll and G HH groups during the first 2- 
year period, controlling for those confounding factors which may 
influence 25(OH)D and muscle mass and strength development 
(age, and the changes of body height, vitamin D intakes, level of 
physical activity and PTH). We further compared the percentage 
change of the muscle mass and strength between the Gll an d 
Ghh groups during the first 2-year period, controlling for those 
confounding factors which may influence 25(OH)D and muscle 
mass and strength development (age, and the changes of body 
height, vitamin D intakes, level of physical activity and PTH). We 
found that, there was an interaction between the menarche status, 
pre- vs. post-groups, and the vitamin D status, sufficient vs. 
insufficient groups, in change of LM WB (p = 0.001), aLM 
(p = 0.003) and mCSA (p = 0.009). By contrast to the Gll, the 
Ghh group increased more in LM WB (p — 0.004, Figure 2a), 
aLM (p = 0.001 Figure 2b) and mCSA (p = 0.027 Figure 2c) 
post-menarche. Within the G H h group, they increased more in 
LMwb (p = 0.005, Figure 2a) and mCSA (p = 0.003 Figure 2c) 
in post-menarche than pre-menarche, while within the Gll group, 
they increased more in LM W b (p = 0.028, Figure 2a) and aLM 
(p = 0.00 1 Figure 2b) during pre-menarche than post-menarche. 
No significant differences between the vitamin D sufficient and 
insufficient groups were found in the percent change of MVC e ib ow 
(Figure 2d) and MY( ;,„., (Figure 2e). However, the G H h group 
tended to increase more their MVC elbow post-menarche compared 
to the Gll group (p = 0.091) and pre-menarche (p = 0.051). 

We found that girls with Apal allele of AA had higher LMwb, 
aLM and mCSA than their counterparts with Aa and aa alleles at 
the baseline, 2-year and 7.5-year follow ups (Table SI). After 
controlling for the level of 25(OH)D the results remained the same. 

Discussion 

Our data demonstrate that, during growth, there was a 
temporary negative association between 25(OH)D and muscle 
mass and strength development prior to menarche due to fast 
growth. After menarche these inverse associations disappeared. At 
the age of 18 years, no differences in muscle mass or strength 
between the 25(OH)D groups were found at the 7.5-year follow- 
up. 

Vitamin D status of children in transition through puberty is an 
important area of study, especially at latitudes above 37 degrees 
north wherein the cutaneous synthesis of vitamin D is only possible 
for no more than half the year [16-18]. The high prevalence of 



25(OH)D insufficiency (<50 nmol/1) in our study is consistent 
with earlier reports showing that 25(OH)D insufficiency is 
common in growing Finnish children and adolescents [19,20]. 
Similarly, a number of studies worldwide, even in sunny areas 
[21], have reported a high prevalence of 25(OH)D insufficiency in 
healthy (asymptomatic) growing children and adolescents [22,28- 
30]. The high prevalence of 25(OH)D insufficiency during 
childhood is alarming and calls for additional research to establish 
a level of 25(OH)D that optimizes skeletal development during 
growth and prevents potential negative health effects later in life. 

Although the role of vitamin D in skeletal health and 
development has been a hot topic, the impact of vitamin D on 
muscle mass and function during growing period is less clear. 
Cross-sectional studies have reported positive relationship between 
25(OH)D, and forearm muscle (grip) strength [31] and jump 
height, velocity and power [32] in adolescent girls. El-Hajj et al. 
[33] showed in a randomized controlled study that vitamin D3 
(cholecalciferol) supplementation significantly increased lean mass 
in 10-17 years old females. Whereas another prospective study 
showed an age-related decline in serum 25(OH)D level in girls 
aged 4-8 years over a period of 1-7 years [32]. The authors found 
that the decline of 25(OH)D level with age was eliminated after 
controlling for LM, suggesting that the increase in LM during 
growth is associated with a decline in circulating 25(OH)D 
concentration. Further, Ward et al. found no improvements in 
muscle force or power with one year vitamin D2 (ergocalciferol) 
supplementation in postmenarchal 12- to 14-yr-old girls [34]. 
Goswami et al. in turn showed that while 6 months cholecalcif- 
erol/calcium supplementation effectively increased and main- 
tained serum 25(OH)D, the increase did not improve skeletal 
muscle strength in young girls [35]. Our data show that the low 
LM WB and aLM, MVC e n,o W an d MVCknee at the baseline in the 
vitamin D sufficient group disappeared at 2-year follow-up 
indicating that the associations between 25(OH)D and muscle 
mass and strength is maturation stage dependent. Moreover, the 
temporary negative association between 25(OH)D) and muscle 
mass and strength prior to pre-menarche does not have long-term 
effects on development of muscle mass and strength as indicated 
by the fact that no differences in muscle mass or strength between 
groups were found when girls were at the age of 18 years. 

Vitamin D receptor (VDR) polymorphisms have been suggested 
to be associated with muscle strength, but evidence for this is 
limited and conflicting [36] . In our study, no significant differences 
in terms of distribution of different VDR genotypes between 
groups were found. However, we found that the girls with VDR 
Apal AA homozygotes of frequent alleles had higher LMwb, aLM 
and mCSA than their counterparts with Aa and aa alleles at 2 and 
7.5-year follow-up, suggesting that certain VDR polymorphisms 
may affect muscle mass development in girls during growing 
period. The molecular basis for the association between VDR 
Apal polymorphism and muscle mass in the present study remains 
unclear, and more research is needed to confirm and explain these 
associations. 

Vitamin D intakes and other lifestyle factors have also been 
shown to influence both level of 25(OH)D [37] and muscle mass 
and strength in children [38] and old women [39] . We found no 
significant differences between the vitamin D insufficient and 
sufficient groups in terms of physical activity and micronutrient 
intake except for the intakes of vitamin D: the insufficient group 
had less consumption of vitamin D than the sufficient group in 
post-menarche. However, the intake of Vitamin D was not 
associated with the level of 25(OH)D and the change of muscle 
mass was independent of change in vitamin D intakes in this study. 
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Figure 2. Comparison of changes of lean mass (a=LM WB & £>=aLM), muscle cross-sectional area (c=mCSA) and muscle strength 
(Gr=MVC e | bow & e=MVC knee ) during the 2-year period in pre- and post-menarche girls at the 2-year follow-up time point. Estimated 
mean with SE (error line) controlled for age and change of body height, vitamin D intakes, level of physical activity and PTH. LM WB = whole body lean 
mass; aLM = appendicular lean mass; mCSA = muscle cross-sectional area; MVC e | bow = maximum voluntary muscle contraction of elbow flexors; 
MVC knee = maximum voluntary muscle contraction of knee extensors; G HH = consistently vitamin D sufficient; and G LL = consistently vitamin D 
insufficient. 

doi:1 0.1 371 /journal.pone.00821 24.g002 



PTH may also have a direct effect on skeletal muscle because 
PTH has been shown to impair energy production, transfer and 
utilization in skeletal muscle in rats [40] and influences skeletal 
muscle protein and amino acid metabolism in rats [41]. PTH is 
also known to increase free intracellular calcium concentrations in 
muscle tissue [42] which may disrupt muscle structure or function 
[43]. However, only a few studies have shown that PTH may 
influence muscle mass and strength, independent of either 
25(OH)D or l,25(OH)2D levels [44,45]. Our data showed that 
high levels of PTH were associated with low levels of 25(OH)D. 
However, the level of PTH was not associated with muscle mass 
and strength gain and the relationship between 25(OH)D and 
muscle mass was independent of PTH. 

Strengths and Limitations of Study 

Despite several strengths of our study, including the longitudinal 
data followed girls from pre-puberty to early adulthood and use of 
hierarchical (multilevel) model in which allowed us to evaluate the 
change of 25(OH)D and muscle mass/ strength during the entire 
growth period. The results of our study must be interpreted in the 
light of its limitations. First, few subjects in our cohort had very low 
levels of 25(OH)D (especially at baseline, 6% <25 nmol/1), 
limiting statistical power to detect an effect of severe vitamin D 
deficiency. Second, tissue levels of vitamin D may be more 
important with respect to muscle mass or strength than the serum 
levels assessed here, although there is now controversy about 
whether skeletal muscle even possesses receptors for vitamin D 
[46]. Thirdly, participants were Finnish girls, and the degrees to 
which genetics, and latitude, or other environmental factors may 
have contributed to these differences remain unclear, hence the 
results may not necessarily apply to other populations. Further- 



more, analyses were adjusted for multiple factors, but the 
possibility of residual confounding factors cannot be excluded. 
Finally, although most of relevant research suggests that circulat- 
ing 25(OH)D is recognized as the best functional indicator of 
vitamin D status, l,25(OH) 2 D was not assessed and this must be 
considered as a limitation. Focusing only on serum 25(OH)D, 
information on the functional role of vitamin D during 
peripubertal muscle growth may not have been captured. 

Conclusion 

Our results showed that when comparing the levels of serum 
25(OH)D between the two groups with lower to higher than 
50 nmol/1, no significant negative influence of 25(OH)D on the 
development of muscle mass and strength during pubertal growth 
were found both with longitudinal and cross-sectional comparison. 
On the contrary, our results suggest that the temporary negative 
association between 25(OH)D and muscle mass arises as a 
consequence of fast growth prior to menarche, and that this 
negative association is diminished through catch-up growth after 
menarche. 

Supporting Information 

Table SI Comparison of 25(OH)D, muscle mass and 
strength among different genotypes in girls at the 
baseline, 2-year and 7-year follow-up (ANOVA by Sidak 
method for multiple comparison; quoted values are 
mean and SD). LMWB = whole body lean mass; aLM = appen- 
dicular lean mass; mCSA = muscle cross-sectional area; MVCel- 
bow = maximum strength of elbow flexors; MVCknee = maximum 
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strength of knee extensors. *<0.05 pairwise comparisons between 

AA and Aa, and AA and aa. 
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